Shell Model Description of Isotope Shifts in Calcium 
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Isotope shifts in the nuclear charge radius of even and odd calcium isotopes are calculated 
within the nuclear shell model. The model space includes all configurations of nucleons in the 
2s, 1^3/21 1/7/21 and 2p3/2 orbits. The shell model describes well the energies of the intruder states 
in Sc and Ca, as well as the energies of the low-lying 2'*' and 3~ states in the even Ca isotopes. 
The characteristic features of the isotope shifts, the parabolic dependence on A and the prominent 
odd-even staggering, are well reproduced by the model. These features are related to the partial 
breakdown of the Z = 20 shell closure caused by promotion, due to the neutron-proton interaction, 
of the ds shell protons into the fp shell. 
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Introduction The appearance of shell gaps associated 
with magic nucleon numbers is one of the cornerstones in 
nuclear structure. However, it has become increasingly 
evident in recent years that these magic numbers, and 
the corresponding shell closures, might get eroded with 
increasing neutron excess. To understand the origin of 
this erosion and to identify whether this trend applies to 
all magic numbers is one of the great challenges in nuclear 
structure with, among other things, strong astrophysical 
implications. A prominent example is the magic neutron 
number = 20 which vanishes in proton-deficient nu- 
clei with Z < 12 . This erosion of the shell closure 
has been related to cross-shell proton-neutron interaction 
which correlates the 2s and 1^3/2 proton orbitals with the 
1/7/2 2p3/2 neutron orbitals and leads to appreciable 
deformation H . If cross-shell correlations are indeed the 
mechanism for the shell erosion, then first indications are 
already visible in the stable calcium isotopes, as we will 
argue in this Letter. Our argument is based on the un- 
derstanding and explanation of the nuclear charge radii, 
(r^), in the calcium isotopes which by itself have been a 
challenge for nuclear theory for the last two decades. 

Experiments based on optical isotope shifts as well as 
on muonic atom data (see Ref. p| for a summary and 
original references, in particular [Ej ) , revealed a char- 
acteristic parabolic shape of the isotope shifts with a 
pronounced odd-even staggering when neutrons fill the 
/7/2 orbit and the mass number changes from A = 40 to 
A = 48. Thus, the challenge consists of understanding (i) 
why ^^Ca is not larger than ^°Ca, {ii) why the maximum 
of {r1) is reached in ''''Ca, {Hi) why odd A calcium iso- 
topes are considerably smaller than the average of their 
even neighbors, and {iv) last but not least, the actual 
magnitude of the isotope shifts. 

There have been numerous attempts to explain these 
findings. If they are indeed related to the cross-shell 



proton-neutron correlations, one expects that in order to 
describe them, one must include higher than 2p-2h corre- 
lations. This is necessary in order to reproduce the odd- 
even staggering and, at the same time, shell effects re- 
lated to the neutron shell closures at = 20 and 28. This 
conjecture is supported by the fact that mean field calcu- 
lations, which usually aim at describing nuclear masses, 
deformation parameters, and radii over a large region of 
nuclear masses and charges, cannot account for the de- 
tails of the calcium isotope shifts. The dependence of 
{r1) on A is usually featureless. Some of the approaches, 
however, are able to account at least for the property (j), 
i.e., the near equality of (r^) in "^"Ca and ''^Ca. To this 
category belong calculations based on the Hartree-Fock 
method with Skyrme interactions j^, relativistic mean 
field methods Q , and the extended Thomas-Fermi model 
with the Strutinsky- integral 

An extension of the mean field method with pairing, 
the so-called local energy- density functional approach , 
can describe most of the features mentioned above by 
introducing three- and four-body forces via the density 
dependence of pairing, and adjusting the corresponding 
free parameters. 

That three-body forces can do the job has been rec- 
ognized long time ago. In |^ the isotope shifts in Ca 
were quite well described using the configuration mixing 
within the (/7/2rf3/2)"' shell model space, and introducing 
an isoscalar particle-hole force whose strength was pro- 
portional to the number of valence nucleons. The origin 
of that dependence and the magnitude of the correspond- 
ing coupling constant, however, remained unexplained. 

Below we show that the modern shell model can de- 
scribe the isotope shifts in calcium consistently and with- 
out the need for many-body forces or adjustable param- 
eters. The basic mechanism of the A dependence of the 
charge radius (r^) remains the same, however, as in Ref. 
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Due to the configuration mixing across the Z = 20 
shell boundary, protons are lifted from the sd to the fp 
shell, resulting in the increase of (r^) equal to 



Srl{A) = ^n%{A)-b- 



(1) 



where Z = 20, 6 is the oscillator parameter which we 
assume remains constant for A = 40 — 48, and the num- 
ber of protons in the fp shell, nJp{A), is the calculated 
quantity. Below we show how it is is calculated and how 
one can relate it to other manifestations of breaching of 
the Z = 20 shell boundary. 

Shell model Progress in the application of modern 
nuclear shell model has been facilitated by the develop- 
ment of numerical codes (the m-code ANTOINE and the 
J-coupled code NATHAN Q ) that can be used with rel- 
atively large valence spaces. The selection of the appro- 
priate single-particle space, and the corresponding effec- 
tive interaction (in particular the choice of its uncertain 
monopole part ) is the most important starting point 
of any shell model application. 

Since here we are interested in the description of cal- 
cium isotopes, it is imperative to include states in the 
vicinity of the N — Z ^ 20 shell boundary. Therefore, 
the chosen valence space consists of the £^3/2, Si/2: /7/2 
and P3/2 subshells for both protons and neutrons. (Thus 
^^Si represents the inert core.) This valence space, first 
used in Ref. has the advantage that the existing 

codes make it possible to describe all Ca isotopes with- 
out truncation. (The largest dimension, 34,274,564 in 
the m-scheme basis, is encountered for the ground state 
of ^^Ca.) The other advantage is the essential absence of 
the spurious center-of-mass motion (no attempt to sup- 
press it has been made). 

For this valence space we start with the two-body ma- 
trix elements (TBME) of Ref. Q , which are defined with 
respect to the ^^O core. The single particle energies are 
now modified to reproduce the ^^Si spectrum. The inter- 
action of Ref. jlj] was built in blocks {sd~sd, sd—fp and 
fp— fp), which incorporate the higher-order excitations 
(core polarization, 2p2h excitations etc.). In particular, 
the 2p2h effects consist mainly in pairing renormaliza- 
tions. As we will include explicitly such mixing, and 
to avoid double counting, a schematic pairing hamilto- 
nian was subtracted from the interaction. Finally, the 
cross sd — fp monopoles of Ref. jlj] were adjusted to the 
masses of neutron-rich isotopic chains. They were here 
retuned to reproduce the gap at ""^Ca and the spectra of 
and 41 Ca. 

Results As indicated in Eq. (1) above, we assume 
that the main cause of variation of the charge radius with 
changing neutron number is the lifting of protons from 
the sd shell into the larger fp shell. This feature then 
suggests that there should be a correlation between the 
isotope shifts and other manifestations of the incomplete 
Z = 20 shell closure, e.g. the appearance of low-lying 



intruder states. In Fig. ^ the experimental excitation 
energies of the intruder 3/2+ states in odd- A Sc isotopes 
are compared with the calculated ones. The agreement is 
very satisfactory, and in particular the lowering of these 
states in ^^Sc and "^^Sc is properly described. 
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FIG. 1. Excitation energies of intruder states in Sc and in 
even Ca isotopes. The experimental (circles) and calculated 
(squares) energies of the 3/2+ states in Sc and the experimen- 
tal (stars) and calculated (crosses) energies of the excited 
states in Ca are shown. 

The excited Oj states in the even Ca isotopes also be- 
long to the intruder category. Their excitation energies 
are shown in Fig. |l| as well. Again, the agreement be- 
tween the experimental and calculated energies is quite 
good. (The larger discrepancy in ^^Ca is presumably 
caused by the absence of the neutron f^ /2 subshcU in the 
model space.) These Oj states are dominantly 4p — 4/i 
states involving both neutron and proton excitations. 
The parabolic behaviour of both intruder excitations in 
Fig. |] suggests that the configuration mixing between 
the naive shell model configurations and the intruders 
involving proton excitations from the sd shell will likely 
be most pronounced near A = 44 where the correspond- 
ing states are lowest. 

The occupation numbers for both protons and neu- 
trons, n/p(A), in the ground states of all A = 40 — 
48 Ca isotopes were then evaluated in the full shell 
model space. It is illustrative to consider, however, how 
these numbers change when more and more nucleons 
are allowed to be excited from the naive configuration 
(ds)^''+^'^, {fp)^^~'^^^" ■ Varying the number of nucleons, 
t, allowed to be excited from the naive configuration 
above, from t = 2 to its maximum value t = A — 2% 
we find that the number of holes in the ds shell satu- 
rates very slowly. For example, in ''^Ca and "^^Ca one 
has to allow t = 6 and 7, respectively, just to reach the 
halfway point of the final n/p(j4), and i = 10 and 11 
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to achieve saturation. This means that the "hfting" is 
a comphcated process, involving substantial rearrange- 
ment of many nucleons. This can be seen also in the fact 
that approximately equal numbers of protons and neu- 
trons are lifted from the ds shell. It is therefore difficult 
to identify a simple cause, or a definite component of the 
hamiltonian, as the driving force of this effect. 

In the full space the quantity nfp{A) is 1.10 in ■*°Ca, 
reaches its maximum of 1.97 in '*"'Ca, and then decreases 
again to 0.78 in "^Ca. 
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FIG. 2. Isotope shifts in calcium. The experimental data 
(circles connected by a solid line) and the shell model results, 
Eq.(l), (stars connected by a dashed line) are shown. 

Knowing the quantities Ufp one can calculate the iso- 
tope shifts from Eq. (1) by subtracting the correspond- 
ing Sr'^{A = 40) from Sr'^{A). (Nucleon form fac- 
tors contribute negligibly to the isotope shift and thus 
are neglected.) The results, calculated with a constant 
b — 1.974 fm as in are shown in Fig. || and com- 
pared with the experimental values. The trends, i.e. the 
properties (i) - (iii) are clearly well reproduced, but the 
magnitude of the calculated shifts is smaller than the ex- 
periment suggests. 

This cannot be attributed to our choice of the oscil- 
lator parameter b. Neither it can be cured by replacing 
the harmonic oscillator wave functions by the Woods- 
Saxon single-particle wave functions, which changes the 
isotope shifts shown in Fig. ^ only very little, provided 
the core (r^) is assumed to be independent of A. (Al- 
though assuming a very slight increase of b{A), i.e. a 
gradual increase of the core (r^) from ^ = 40 to ^ = 48, 
might "straighten" the curve a bit.) Moreover, choos- 
ing the recommended average dependence of b on the 
mass number A and isospin T ||l|] (6^ = 1.07^^/^(1 - 
{2T/Af)exp{3.5/A) fm^) would make the radius of ^^^Ca 
considerably larger than the radius of ^°Ca. The b{A, T) 
dependence, which on average follows Q, clearly is sub- 



ject to local shell effects. Most likely, the calculated am- 
plitude of the isotope shift might be modified if the so far 
neglected effects of the filled and empty orbits 
are added. 

Collective states When considering isotope shifts, one 
also has to take into account the effect of the zero-point 
motion associated with the surface vibrational modes (see 
p5[ for a systematic approach to this question. There 
it is stressed that only the low-lying low-multipolarity 
surface vibrational states contribute significantly). The 
collective states increase the mean square radius by 



A 



(2) 



where i?o is the equivalent sharp nuclear radius, and P\ 
is the vibrational amplitude of the mode A, related to the 
corresponding B{EX,0 J = \) value. What matters 
for the isotope shifts, naturally, is not the absolute value 
of the amplitudes (3\, but their variation with A. 




FIG. 3. Excitation energies of the 2"*" and 3 states in the 
even Ca isotopes. 

Formula was derived as a correction to the mean 
square radius calculated using the mean-field methods. 
On the other hand, since the shell model includes all 
correlations of nucleons in the valence shells, blind use of 
it would amount to at least partial double counting. 

The energies of the lowest lying 2+ and 3^ states are 
compared to the shell model in Fig. |^. The good agree- 
ment there shows not only that the chosen interaction de- 
scribes this important spectroscopic data well, but that 
the zero-point motion effect on the nuclear radius is (at 
least most of it) automatically included. 

In addition, in Fig. ^ we compare the experimental and 
shell model B{E2, 2-^0) values. (The calculations were 
performed with the usual effective charges Cp = 1.5, e„ = 
0.5.) Again, the agreement is quite good, strengthening 
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FIG. 4. Comparison of the experimental and shell model 
B{E2, 2-^0) values in the even Ca isotopes. 



our belief that most of the effect of the collective states 
is already contained in our calculation. 

Discussion The present calculation shows that the 
Z — 20 and N — 20 shell boundaries are clearly not 
absolute, and that substantial configuration mixing in- 
volving nucleons from the ds shell is present even in the 
Ca isotopes. 

This finding is in line with the previously recognized 
and perhaps more dramatic "islands of inversion" related 
to the configuration mixing involving the = 20 shell 

Here 



boundary in neutron rich nuclei (see Ref. |16 



we are encountering a similar situation, now however in 
stable nuclei near or at iV = Z. 

Is there any other evidence for the incomplete shell clo- 
sure in '^''Ca? Indeed, the study of the {n,p) and {p,n) 



reactions ]1^20| on Ca revealed that the integrated 
Gamow-Tcller strength B{GT) = 1.6±0.1 below 15 MeV. 
Obviously, if the ds shell in '^'^Ca were really closed, the 
B{GT) would vanish. Note, that without the full in- 
clusion of the /5/2 and ^5/2 orbit we can only estimate 
the total Gamow- Teller strength in ^''Ca. Without these 
orbits we obtain the total i3(GT)=0.53 using the usual 
quenching of 0.744 |2^. When we allow the t = 1 (i.e. 
one particle or one hole) excitations involving the 
and (i5/2 orbits, we obtain B{GT) = 2.75 with quench- 
ing, in fair agreement with the experiment. Clearly, some 
of the strength associated with these orbits will be above 
15 MeV, and should not be counted. 

In conclusion, the present calculation shows that the 
isotope shifts, and the position of the intruder states in 
Ca isotopes, can be well described by large scale shell 
model calculations. The judicious choice of the valence 
space and the corresponding effective hamiltonian is the 
key ingredient to this success. Thus, the challenge to the 
nuclear structure theory, described in the introduction 



to this work, has been largely met. No new forces or 
parameters are needed to describe the dependence of the 
nuclear radius in the odd and even Ca isotopes on the 
mass number A. 
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